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Abstract 

Spatiotemporal chaos arising from the competition between sine-Gordon- 
breather and kink-antikink-pair sohtons by reshaping an ac force is demon- 
strated. After introducing soliton collective coordinates, Melnikov's method is 
applied to the resulting effective equation of motion to estimate the parameter- 
space regions of the ac force where homoclinic bifurcations are induced. The 
analysis reveals that the chaos-order threshold exhibits sensitivity to small 
changes in the force shape. Computer simulations of the sine-Gordon system 
show good agreement with these theoretical predictions. 

I. INTRODUCTION 

Although nonlinearity is an ubiquitous and fascinating feature of real-world dynamical 
phenomena, some aspects of such a fundamental property remain poorly explored, partly 
because they are rarely incorporated into the mathematical models which aim to describe 
even the simplest of such phenomena. In the context of nonautonomous systems, one of these 
aspects is the nature of the periodic excitations which can act on a given nonlinear, dissipa- 
tive or Hamiltonian, system. Since nonlinear systems represent the general case, it would be 
natural to model periodic excitations (external or parametric) by using periodic functions 
that are solutions of (physically relevant) nonlinear equations, instead of harmonic functions 
which are solutions of linear systems, even though the latter have been overwhelmingly used 
up to now. To consider nonlinear periodic excitations implies enlarging the amplitude-period 
parameter space to include the parameters, ctj, that control the excitation shape. Note that, 
in the general case, these parameters a-i are the (normalized) Fourier coefficients of the pe- 
riodic function. In physical terms this means that, for a fixed amplitude and period, the 
parameters Ui are responsible for the temporal rate at which energy is transferred from the 
excitation mechanism to the system. In the context of systems described by ordinary dif- 
ferential equations, this idea has led to the demonstration of the existence of new generic 



temporal chaos*->order routes by changing only the shape of a nonhnear periodic excitation 
[1-3]. 

In this present work, the onset of spatio-temporal chaos associated with the competition 
between sine-Gordon-breathers and kink-antikink-pairs is investigated in a one-dimensional 
sine-Gordon (sG) system 

$tt - $xx + sin $ = -(5$i + 7 sn [4:K{m)t/T; m] , (1) 

when only the excitation shape is varied. Here, sn {u; m) is the Jacobian elliptic func- 
tion (JEF) of parameter m {K{m) is the complete elliptic integral of the first kind) used 
solely as an example to illustrate the reshaping-induced order-chaos route. When m = 0, 
then sn [4K{m = 0)t/T;m = 0] = sin(27rt/T), i.e., one recovers the previously studied case 
of harmonic excitation [4]. This is relevant in comparing the structural stability of the 
sG system when solely the excitation shape is varied. In the other limit, m = 1, one 
obtains the Fourier expansion of the square wave function: sn[4K{m = l)t/T;m = 1] = 
(4/7r) Yl'^=o (2^ + ^y ^^^ [(2^ + ^) Svrt/T]. From an experimental standpoint, the choice of 
the above JEF is motivated for the fact that only the excitation shape is varied by increasing 
m from to 1, and there is thus a smooth transition from a sine function to a square wave. 
While these two limiting shapes are readily implemented in experiments, no theoretical 
analysis connecting them has been developed up to now in the context of sG systems. 

The organization of the paper is as follows. In Section 2, a collective coordinate approach 
[4] is used to obtain an effective equation of motion for the collective variables. In Section 
3, Melnikov's method is applied to the resulting effective equation of motion to estimate 
the parameter-space regions of the elliptic ac force where the onset of chaotic instabilities is 
expected. Numerical evidence supporting the theoretical predictions from previous sections 
is presented in Section 4. Finally, Section 5 includes a summary of the findings. 



II. COLLECTIVE COORDINATES 

In this present work, a single-breather excitation is assumed to be present in the absence 
of dissipation and forcing (5 = 7 = 0). It is well known that the simultaneous action 
of dissipative {6 > 0) and spatially uniform harmonic forces (7sin (27rt/T)) can break the 
breather into a kink-antikink pair, which can then recombine into a breather soliton [4,5]. For 
fixed dissipation, this process may occur repeatedly yielding an intermittent chaotic sequence 
for certain values of amplitude and period. The above elliptic ac driving permits one to study 
such a competition between the breather and the kink-antikink pair for parameter values 
belonging to /arg'er regions of the amplitude-period parameter space. To this end, one makes 
a severe mode truncation to the breather collective coordinates [6,7], which yields a set of 
ordinary differential equations capable of being studied with the aid of Melnikov's method 
[8,9]. Such an effective equation exhibits a separatrix dynamics which is equivalent to the 
separatrix dynamics of the full sine-Gordon equation. In particular, homoclinic bifurcations 
correspond to (chaotic) intermittent binding and unbinding of the kink-antikink pair. In 
order to determine the equations of motion for the breather collective modes, one first 
needs the breather and kink-antikink solutions to the unperturbed sine-Gordon equation 
(5 = 7 = in Eq. (1)) [4]. 

Breather: 



1 f tan ^ sin \v (t)] 1 



v(t) = vq + tcos9, 
k = - sin 9, 
E = 16 sin 9. (2) 

Kink-antikink pair {K K): 



Uj^f-.{x,t) = 4 tan 



-1 



sinh[M(t- to) /(!-«')'/'] 



^^'■^' ■'^ ^ ^"^^ 1 U cosh [{X - Xo) /(I - m2)1/2] 

16 

u 



^i^x-TT— ;iW^>16,«^0. (3) 



Here, u is the t ^ oo velocity of the kink while —u is the t ^ oo velocity of the antikink, 
and vq, Xq, 6 {0 ^ 6 ^ 27r) are constants. Also, ujb = cos^ is the breather's internal frequency 
and E, E^^^ are the breather and kink-antikink-pair energies, respectively. Following Ref. 
[4], one assumes the ansatz that the solution of the perturbed sine-Gordon equation has 
the same form as the unperturbed breather (2), but now allowing 9 and vq to be temporal 
functions: 

v{t) =VQ{t)+ f cos 9 {t')dt'. (4) 

This approximation means that one is assuming a perturbation sufficiently weak so that its 
principal effect on the system is to continuously alter the frequency and the phase of the 
breather. Therefore, such an approximation is valid when ujb ^ 1 since the breather energy 
is Eb = 16(1 — lUbY^'^- Next, one can deduce the equations of motion for the breather 
collective modes {u,v), where 



u(t) = Ub{x = Xo,t) = 4tan"M(t), 



OUe 
dt 



v{t) 

L ^^ Jx=xo 

A{t) = tan 6* sin t; (t) . 



4:dA{t)/dt 
l + A{tf' 



(5) 



After the coordinate transformation z = tan ('u/4) , z= w, such equations of motion can be 
shown to be 



X =f(x) + eg(x,V'), 
ij =AK{m)/T, 



(6) 



with 
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f(x) 
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-zfl-^i;2^Q^^2^-l 



(7) 



and the perturbation 

/ 

eg (x; t) = 
V 

where 





(8) 

h{t)gp{z,w) -5gsiz,w) 



h{t) = -'jsn[AK{m)t/T;m\ , 
, , Uw^ + z^)z^-s{z^ + w^z-' + {l-w^){l + zm] 

q [2^ -y; J ^ i-^ 1: i __ 

{l + sw^){w^ + z^){l + s){l + z^f'^ 
2w^z^s^ (1 + z^) + w^s^ (1 + 2z^) (1 + z^)' 

(1 + SW2) (^2 + ^2) (1 + 5) (1 + ^2)1/2 ' 

w [(2 + ^=^) (w^ + z^) z^ -s{z^- w^ [(2 + z^) z*-{l- w^) (1 + z'')]}] 
^^^^''"^~ {1 + sw^) {w^ + z^) {1 + z^) 

w [wh^ (1 + 2^2) (1 + Z^)] 



{1 + sw"^) {w^ + z^) {1 + z^)' 
^ smh-\z) 

zil + z^f'- ^' 

Noting that the separatrix in the phase-plane u — v is given by 

Mo,±(t) = ±4tan"^(t), 

t;o,±(t) = ±4/(l+t2), (10) 

one finds that on the separatrix in the phase-plane z — w 

Zo,i{t) = ±t, (11) 

^i^O,±W = ±1, 

where the positive (negative) sign refers to the upper (lower) homoclinic orbit in the phase- 
plane u — V, which is obtained for the energy i? = 16 of the unperturbed energy functional 
E = 2tan\u/4) + {u /A^ sec'^{u/4). 

III. CHAOTIC THRESHOLD FUNCTION 

Now one can estimate the threshold of chaotic instabilities for the effective system (6) 
using Melnikov's method. Thus, the corresponding Melnikov function is 
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M±(to) = tD + CY^ an{m)pn{T) sin 



n=0 

where 



(2n + 1) 27rto 
T 



(12) 



oo 



fl — tanh^ x] \{2 — tanh^ x) tanh^ a; + r + r^ f 1 — tanh^ x)l 
Z^ = 5 / -!^ '-^ '- 5^ 'Mx ~ 2.03727(5, 

' 1 + r 



— oo 

oo 



COS flS^ sinha;) Tfl — r) fl — tanh^x)! sechx + 2r^ tanh^a; + r^ fl + tanh^x) 

p (2,/r) . j^ -^ U^^^ JL^ ^ i... 

C = 7r27/(4v^is:(m)), 
a„(m) = csch [{2n + 1) TrfsTfl - m)/{2K{m))] , 
p^iT)=p[{2n+l)27r/T], (13) 

and where r = 2a: csch (2a;). From Eq. (12) one sees that a homochnic bifurcation is 
guaranteed if 

-<f/(m,T), (14) 

7 

where the chaotic threshold function is 

^K^)^ , /^^, J ^(-l)"aJm)p„(T). (15) 

^ ^ ' n=0 

The series in Eq. (15) satisfies the conditions of the Leibnitz theorem on alternate series, and 
then one can deduce a slightly more restrictive condition for the onset of the aforementioned 
spatio-temporal chaotic instabilities, in order to facilitate the comparison with the previously 
studied case [4] of a harmonic force {m = 0): 

r 

-<U{m = 0,T)R{m), (16) 

7 

where U{m = 0,T) = 7rp(27r/T)/(4a), i?(m) = tt {^K{m)y^ csch [nKil - m)/{2K{m))]. 
It is worth noting that the condition (16) is a reliable approximation due to the coefficients 
a„(m) ~ 0, Vn ^ 1 and Vm G [0,771*] with m* being very close to 1. A plot of the shape 
function R{m) is shown in Fig. 1, where one sees that R{m = 0) = 1, as expected. The 

7 



rapid growth of R{m) as m -^ 1 means that the chaotic threshold exhibits hypersensitivity to 
changes in the excitation shape as this approximates a square wave. This can be understood 
in terms of the increasing accumulation of effective harmonics in the Fourier series of the 
elliptic force as m ^1. Note that this number of effective harmonics is not a linear function 
of m because of the dependence of the elliptic force on K{rn). 

IV. NUMERICAL SIMULATIONS 

For the purpose of the comparison with the above theoretical predictions, the parameters 
5 = 0.2, T = 27r/0.6,and L = 24 were fixed, and an sG breather with frequency cub = v2/2 
and Xo = f = (cf. Eq. (2)) was chosen as initial data. The period T = 27r/0.6 corresponds 
to a frequency uj = 0.6 (for a sinusoidal force, m = 0) for which the chaotic attractors are 
dominated by breather to kink-antikink transitions [11,12]. Equation (1) was integrated 
under periodic boundary conditions $ (a; = —L/2,t) = $ (a; = L/2,t) for all t. Figure 2 
shows an illustrative example of the comparison between the theoretically obtained chaotic 
threshold amplitude '-fth (solid line, cf. Eq. (16)) and the corresponding numerically obtained 
amplitude (dots) as the shape parameter m varies. One sees that the amplitudes obtained 
numerically are slightly larger than the respective theoretically predicted ones, as expected, 
although an overall good agreement is found. Similar good agreement was found for other 
sets of parameters in the large-period regime (T ^ 27r), provided that the damping and 
driving force were sufficiently weak. An illustrative sequence of the order-chaos route is 
shown in Fig. 3 where the initial state (m = 0) occurs after the system has undergone a 
spatial period doubling (for m = and 0.6 ^ 7 ^ 0.9, cf. Ref. [12]) but remains simply 
periodic in time. This two-breather state has thus a spatial period | and a temporal period 
1, and exhibits mirror symmetry with respect to the origin. The spatial symmetry is broken 
as the forcing shape deviates from sinusoidal. FigureS illustrates this phenomenon with 
the values m = 0.8 and m = 0.82. This spatial symmetry breaking of the two-breather 
state appears to be an unavoidable precursor of the onset of spatio-temporal chaos via the 



kink-antikink collision mechanism, as can be appreciated in Fig. 3 for the shape parameter 
m = 0.898 ~ iTLth, which corresponds to the theoretically predicted chaotic threshold (cf. 
Eq. (16)). For higher values of the shape parameter, i.e., when the driving shape is ever 
closer to a square wave, one finds the persistence of the chaotic dynamics (cf. Figs. 3, 
m = 0.99) predicted by the Melnikov-method-based analysis. 

V. CONCLUSIONS 

In sum, it has been shown through the example of an elliptic ac driving that the competi- 
tion between sine-Gordon-breather and kink-antikink-pair solitons strongly depends on the 
shape of the spatially uniform ac driving. In particular, it was theoretically predicted and 
numerically confirmed that a square wave is more effective than a sinusoidal function to yield 
spatio-temporal chaos associated with the competition between sine-Gordon-breathers and 
kink-antikink-pairs. Potential application of the present results to reshaping-induced con- 
trol (suppression and enhancement) of chaotic dynamics [1,14,15] range from DNA molecules 
[16] to long Josephson junctions and quasi-one-dimensional ferromagnets [17]. An interest- 
ing problem under current investigation is whether the above predictions on the chaotic 
threshold remain reliable in the so-called (for m = 0) nonlinear (cubic) Schrodinger regime 
(i.e., T < 27r), where the chaotic dynamics is dominated by breather-breather transitions 
instead of the breather to kink-antikink transitions of the regime studied here (T ^ 27r), 
since both mechanisms involve homoclinic crossings [13]. 
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1. Figure Captions 

Figure 1. Shape function R{m) (see the text) versus the shape parameter m. 

Figure 2. Theoretical (solid line, cf. Eq. (16)) and numerically obtained (dots) chaotic 
threshold amplitude, •jth, versus shape parameter, m, for 6 = 0.2 and T = 27r/0.6. 

Figure 3. Space-time evolution of $ (x, t) for the sG system through four driving periods 
for 6 = 0.2,7 = 0.79, T = 27r/0.6, and shape parameters m = 0, 0.5, 0.8, 0.82, 0.898 ~ rrith, 
and 0.99. 
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